Total cross sections for neutron-nucleus scattering 
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Systematics of neutron scattering cross sections on various materials for neutron energies up to 
several hundred MeV are important for ADSS applications. Ramsauer model is well known and 
widely applied to understand systematics of neutron nucleus total cross sections. In this work, we 
examined the role of nuclear effective radius parameter (ro) on Ramsauer model fits of neutron total 
cross sections. We performed Ramsauer model global analysis of the experimental neutron total 
cross sections reported by W. P. Abfalterer, F. B. Bateman, et. al.„ from 20MeV to 550MeV for 
nuclei ranging from Be to U . The global fit functions which can fit total cross section data over 
periodic table are provided along with the required global set of parameters. The global fits predict 
within ±8% deviation to data, showing the scope for improvement. R has been observed that a 
finer adjustment of ro parameter alone can give very good Ramsauer model description of neutron 
total scattering data within ±4% deviation. The required ro values for Ramsauer model fits are 
shown as a function of nuclear mass number and an empirical formula is suggested for ro values as 
a function of mass number. In optical model approach for neutron scattering, we have modified the 
real part of Koning-Deleroche potentails to fit the neutron total cross sections using SCAT2 code. 
The modified potentails have a different energy dependence beyond 200MeV of neutron energy and 
fit the total cross sections from Al to Pb. 
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The concept of an accelerator driven sub-critical (ADS) 
system is drawing worldwide attention yj, |2j. In this 
ADS system, neutrons are produced by bombarding a 
heavy element target with a high energy proton beam 
of typically above l.OGeV with a current of > 10mA 
[H . Such a system serves a dual purpose of energy mul- 
tiplication and waste incineration. In this context it is 
important to study the systematics of neutron scatter- 
ing cross sections on various nuclei for neutron energies 
up to several hundred MeV. In this paper we propose an 
empirical approach to reproduce the energy and atomic 
mass dependence of total cross sections [a tot ) ■ The nu- 
clear Ramsauer model was first proposed by Lawson [3] 
in the year 1953 as a simple means to understand the 
energy dependence of total cross sections of neutron nu- 
cleus scattering. In order to appreciate this model, it is 
necessary to discuss briefly the optical model (OM) de- 
scription of neutron scattering. In the OM approach, 
the scattering amplitudes are obtained by solving the 
Schrodinger's equation with complex optical model po- 
tentials (OMP), such as in SCAT2 optical model code Q 
with Koning and Deleroche potentials Q. The absorp- 
tion arising from imaginary part of OMP gives the reac- 
tion cross section. The calculations are usually performed 
using partial wave expansion method and the phase shifts 
St ( rji—e 2tS '-—aie l ^ e ) are determined. For a given set of 
potentials, these complex phase sifts are strongly angu- 



lar momentum and energy dependent, 
phase shifts and the wave number (A 



In terms of the 
= H/V2mE) and 



neglecting spins, various cross sections are given below. 
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In Figure[TJ we show the transmission as a function of an- 
gular momentum t for n+ 208 Pb system at high and low 
energies. The results in this figure are generated using 
SCAT2 code |4J with Koning and Deleroche potentials [5j . 
The transmission, which is related to absolute value of 
the complex phase shift (1 — l^l 2 ), is approximately con- 
stant up to some £ m ax for each energy and falls sharply 
after that cut off angular momentum.. We have plot- 
ted the SCAT2 transmissions and seen that this angular 
momentum dependence is very similar in all the systems 
considered in this work. These results are well known in 
the field of heavy ion fusion around the Coulomb barrier 
energies where, a sharp cut off or diffuse cut-off models 
are often used. When multiplied by factor, The transmis- 
sion functions give the well known triangular distribution 
shape of partial wave cross sections ( ai — {it / k 2 ) (21+1) . 
Ramsauer model in the following discussion is essentially 
based on this behaviour of transmission function. 
An excellent experimental data base of neutron total 
cross sections is presently available in the energy range up 



to 600 MeV [6j-|8[ . The most recent work by Koning and 
Delaroche Q presents a very exhaustive search for OMP 
parameters that fit the data very well up to 200 MeV. 
We made a phenomenological Ramsauer model analy- 
sis of the experimental total cross sections. The nuclear 
Ramsauer model [3] provides a simple method to param- 
eterize the energy dependence of neutron nucleus total 
scattering cross sections. This model assumes that the 
scattering phase shifts are independent of angular mo- 
mentum (£) as reflected in Eq.© (rj = ae 1 ' 3 ), in contrast 
to the predictions of the optical model given in Eq.©. 
Further, it was proposed that the ^-independent phase 
shift varies slowly with energy. Peterson applied this 
model for the study of neutron scattering from various 
nuclei [§, LL0| . There were some attempts |lll - fl4| (see ref- 
erences therein) to put this Ramsauer model on a sound 
theoretical basis. The neutron total cross sections have 
thus been well studied |13l - [l9j using this model, over a 
wide range of nuclear masses and neutron energies up to 
500 MeV. In |2CJ,[2l|], simple functional forms for the total 
cross sections through the £q parameter corresponding to 
partial wave distributions were shown. In their method, 
one can determine the partial wave cross sections, con- 
sidering spin. 
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I. ANALYSIS OF NEUTRON TOTAL CROSS 
SECTIONS 

We performed the Ramsauer model fits to the experi- 
mental neutron total cross sections using Eq.©. The 
quantities R(fm), a and j3 are functions of atomic mass 
number (A) and the center of mass energy (E). 
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Figures[2ja,b) show the experimental data represented by 
symbols and the Ramsauer model fits (solid lines) for atot 
cross sections using Eqs. ©- ©. The fits are obtained 
with twelve optimised parameters as given in Eq.© over 
wide mass range of 24 Mg to 208 Pb (Fig. Ufa)) and 9 Be 
to 24 Mg (Fig. 12(b)) (Hf. These fits cover the neutron 
energy region (Ei a b) of 20-550 MeV. Our functional de- 
pendence on energy and mass given in Eq.©-© with 



FIG. 1: Transmission function as a function of angular mo- 
mentum for n+ Pb system, showing the sharp cut off be- 
haviour beyond some £ m ax- The results are generated by 
SCAT2 code using KD potentials. The neutron incident en- 
ergies are indicated on curves. Similar dependence is seen for 
all nuclei considered in Fig[5] 



twelve global parameters ( Eq.©) reproduced the ex- 
perimental data well [24]]. The percentage deviations of 
these fits from data shown in Figs. [Ha,b) have been cal- 
culated. We observed that the fits are within ±6% for 
heavy nuclei and within ±8% for light nuclei for energies 
between 20MeV-550MeV. These deviation functions are 
shown for all systems in FigsJ3(a,b). As seen in Figures 
[3l overall fits are within 8% error from the data. Note 
that the radius parameters for light systems are different 
as mentioned on the Figure [2(b). In order to test these 
parameterization, we predicted a to t using Eqs. (HI)-© for 
six nuclei which were not included in global fits of Fig- 
ure [2J These predictions and the respective experimental 
data are shown in Figure 2J We have seen that the pre- 
dictions are within 8% deviation from data. Model fits to 
total cross sections were already shown by various groups 
[lj, LL31 [20, [21| (see the references therein). Some works 
reproduced very well the Ramsauer peak, superposed on 
their smoothly varying parts of the cross sections [21j . 
The Ramsauer model global fits of experimental neutron 
total cross sections are shown in figures [2J and tested for 
systems which are not used in fitting process as shown in 
FiglU The fit functions and global set of parameters are 
provided in Eqs.©-©, ©. These functions and param- 
eters are a convenient starting point for Ramsauer model 
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FIG. 2: Ramsauer model fits (solid lines) to the experimental 
data (symbols) of total cross sections versus Ei a b using Eq.(|4|) 
[241 ] . The twelve parameters optimised and the functional 
forms used are mentioned in the text in Eqs.([5])-([8]). Fig. (a) 
is for heavy nuclei and Fig. (b) is for light nuclei. The radius 
parameters used for light systems are mentioned on the Figure 
(b). 



analysis of the experimental data and further improved 
fits can be achieved by local finer variation of these pa- 
rameters. In this context, we observed that by adjusting 
only a single parameter tq, it is possible to achieve better 
fits to experimental data, as has been earlier shown by 
us for neutron total scattering on targets 238 U and 232 Th 
[25j ]. We extended this study for all the systems, light and 
heavy considered in figures 121 II It is well known that the 
nuclear radius is given by rgA 1 ' 3 . For Ramsauer model 



FIG. 3: The percentage deviations of the empirical fits in Figs. 
??(a,b) from the experimental data are shown for heavy and 
light nuclei respectively in Figs.(a,b). 



analysis, this is a good approximation only for global 
fits. However, local minor deviations may occur for nu- 
clear radii from this formula roA 1 ' 3 . It is important to 
note that this radius, roA 1 / 3 , should not be taken to 
mean the actual radial size of the nucleus. It represents 
the energy independent part of the radial size required in 
the Ramsauer model fits, see Eq.©. Therefore, it rep- 
resents the size required by neutron total cross sections, 
very similar to the strong absorption radius in heavy ion 
elastic scattering. Hence, this radial size is usually much 
larger than nuclear size. 

Therefore we performed an elaborate study by varying 
the rg parameter and fixing all other parameters as given 
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FIG. 4: The Ramsauer model predictions (solid lines) to the 
experimental data (symbols) of total cross sections versus E; a b 
using Eq.Q, using the functions and parameters that were 
optimised for Figure. ??(a) as in Eqs.([5])-([8|). 











FIG. 6: Ramsauer model fits (solid lines) to the experimental 
data (symbols) of total cross sections versus E; a t using Eq.Q. 
The global fits are same as in Fig[2] The red curve shows 
the fits by varying the ro parameter. Left panel shows the 
fits to cross sections for 238 U and 232 Th systems and right 
panel shows the percentage deviation of respective fits from 
the data. The red symbols show less deviation compared to 
black symbols (global fits) as expected. 
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FIG. 5: The ro parameter values that result in best Ramsauer 
model description to experimental data have been shown ver- 
sus mass number. 




in EqlU We obtained r adjusted fits and their percent- 
age deviations for all the systems considered in this study 
and compared these with the global fits mentioned ear- 
lier. The required ro values are shown in FigjS] and will 
be discussed later. These detailed fits and comparisons 
are shown in Figsl6l- IT6l In these figures, we show cross 
section fits in left panel and percentage deviation of fit 
from data in right panels. In these left panel figures, 
symbols represent data, black curves for global fits and 
red curves for ro adjusted fits. In the right panels, we 
show %deviations of global fits (black circles) and ro ad- 
justed fits (red circles). As can be seen in all these figures, 
the systematic trend in % deviations were removed by ro 
adjustment. The % deviations are in most cases within 
±4%. It is owing to the requirement of a finer adjustment 
of ro parameter values, the global fits produced large de- 
viations up to ±8%. 

As discussed earlier, the ro parameter values that result 
in best Ramsauer model description to experimental data 
are shown in Figj5] The constant value of ro = 1.42988 
of global fits is shown in the figure by black straight line 
for reference. The smooth red solid curve in Figj5] shows 
the required variation of mass dependence of the ro- This 
is given by a formula r (A) = 1.0005 + 0.38045.4 002581 . 
This functional form and parameters show that ro^ 1 ' 3 
is an approximation, however, the finer deviation of nu- 
clear mass dependence is shown by second term. The 
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FIG. 7: As shown in Figs® for 208 Pb and 209 Bi. 



FIG. 9: As shown in Figs© for 93 Nb and n<2t Sn. 
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FIG. 8: As shown in Figs® for 181 Ta and 197 Au. 



FIG. 10: As shown in Figs® for 89 Y and 90 Zr. 



fine adjusted rg values using Ramsauer model are shown 
in Table. 1. in column 2 and the from the formula given 
in column 3. These results in Table are to be compared 
with global fixed value of ro=l- 42990 . Using this depen- 
dence of ro values would result in improved Ramsauer 
model fits of the experimental data. 



II. OPTICAL POTENTIALS FOR NEUTRON 
TOTAL CROSS SECTIONS 

As mentioned in the introduction, the SCAT2 code 
[4[ with Koning Deleroche (KD) [fj} potentials, explains 
consistently all the measured quantities of neutron nu- 
cleus cross sections, such as total cross sections, angular 
distributions, shape elastic cross sections etc. However, 
as well known its use is limited to IKeV to 200 MeV 
energy region only. After achieving good fits using Ram- 
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FIG. 11: As shown in FigsH for 55 Mn and 59 Co. 



FIG. 13: As shown in FigsH for na *Mg and 27 A1. 
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FIG. 12: As shown in Figs(6] for 31 P and 40 Ca. 



FIG. 14: As shown in Figs[6] for 13 C and 19 F. 



sauer model, we have tried to modify the KD potentials 
in SCAT2 program such that they explain the total cross 
sections up to 550 MeV. We use KD potentials up to 
energy E c and a modification function f(E) beyond E c 
for real potentials. Effectively, f(E)=l below E c and the 
f(E) function normalizes the volume term of real part of 
optical potentials, i.e., pote(l) array in SCAT2 code. We 
have chosen exponential energy dependence factor given 
by f(E) = exp((E — E )/a e ), with E and ae and E c be- 
ing adjusted for various nuclei and these values are shown 
in Table [TTJ In the following figures, we show these de- 



tails with and without modification of KD potentials. In 
figures fT7ll20l symbols show the experimental total cross 
sections of W. P. Abfalterer, F. B. Bateman, et. al.,]£$. 
The black lines are for normal SCAT2 code results with 
the prescribed KD potentials with option 7 = 0. This 
option refers to no relativistic corrections in potentials. 
THe blue curves represent results with KD potentials 
with option 7 = 1. As well known, this option refers 
to relativistic corrections implemented to optical poten- 
tials. The red curves show the total cross sections from 
our modified KD potentials using option 7 = 1. As seen 
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FIG. 15: As shown in Figs O for 
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FIG. 16: As shown in FigsEJ for 9 Be. 



in figures, the red curves simulate the neutron total cross 
sections well up to 550 MEV. 



TABLE I: Table shows ro values for various nuclei required 
to achieve best fits using Ramsauer model for neutron total 
cross sections. First column shows nucleus, second column 
shows the ro value. This is to be compared with a fixed value 
of 1.43fm used in global analyses. 



Nucleus 


r (/m) 


formula 


238 tt 


1.4442 


1.43866 


232 pi 


1.44853 


1.43838 


209 B j 


1.43457 


1.43720 


208p b 


1.42689 


1.43714 


197 Au 


1.43534 


1.43558 


181 Ta 


1.43204 


1.43653 


120 Sn 


1.41503 


1.43135 


93 Nb 


1.43409 


1.42817 


90 Zr 


1.4285 


1.42780 


89y 


1.43287 


1.42768 


59 Co 


1.43109 


1.42317 


55 Mn 


1.41083 


1.42241 


40 Ca 


1.43506 


1.41895 


31p 


1.43062 


1.41621 


27 A1 


1.39221 


1.41473 


24 Mg 


1.42302 


1.41347 


19 p 


1.42892 


1.41099 


13 C 


1.38432 


1.40699 


"B 


1.36467 


1.40524 


io B 


1.44432 


1.40425 


9 Be 


1.40193 


1.40315 



TABLE II: Table shows best adjusted -Eo, ae values of the f(E) 
function that fit the neutron total cross sections for various 
nuclei using SCAT2 program. First column shows nucleus, 
second and third columns show the E c ,Eq and a e values. 



Nucleus 


E c (MeV) 


Eo(MeV) 


a e (MeV) 


2U9 Bi 


350 


350 


260 


208p b 


325 


325 


275 


197 Au 


350 


350 


245 


181 Ta 


300 


300 


280 


120 Sn 


250 


250 


270 


93 Nb 


225 


225 


280 


90 Zr 


225 


225 


280 


89y 


225 


225 


285 


59 Co 


225 


225 


275 


55 Mn 


210 


210 


280 


40 Ca 


225 


225 


275 


31p 


250 


250 


270 


27 A1 


225 


225 


285 


24 Mg 


250 


250 


275 


19p 


250 


250 


275 


13 C 


250 


250 


300 



III. CONCLUSION 

In conclusion, we performed the Ramsauer model param- 
eterization of experimental neutron total cross sections 
and derived the systematics of the Ramsauer model pa- 
rameters. The percentage deviation of fits from data were 
plotted and displayed to be within ±8%. It was observed 



that the Ramsauer model fits are sensitive to nuclear ra- 
dius parameter. By fine adjustment of ro parameter lo- 
cally one can obtain very good Ramsauer model fits to 
neutron total cross sections. The r adjusted fits have 
been compared with the global fits for various systems. 
These percentage deviations have been presented in de- 



the total cross sections well up to 550 MeV. 
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FIG. 17: Total cross sections for n scattering on 
13 C, 19 F, 24 Mg, 27 Al. Symbols show experimental data. Black 
curve is SCAT2 results with KD potentaials with 7 = op- 
tion. The blue curve represents with 7 = 1 option and red 
curves represent modified KD potentials together with 7 = 1. 



FIG. 19: Total cross sections for n scattering on 
89 Y, 90 Zr, 93 Nb, 120 Sn. The curves are as explained in previ- 
ous figure. 
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FIG. 18: Total cross sections for n scattering on 
31 P, 40 Ca, 55 Mn, 59 Co. The curves are as explained in previ- 
ous figure. 



FIG. 20: Total cross sections for n scattering on 
181 Ta, 197 Au, 208 Pb, 209 Bi. The curves are as explained in pre- 
vious figure. 



tail. The r parameter values that fit best have been 
shown as a function of nuclear mass. We modified the KD 
potentails incorporating new energy dependence. With 
these new potentials, the SCAT2 code was able to predict 
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